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■iDtCDUCTJCN' 

Tho  continuing  need  for  improved  detection  porfori'ianco  by  surfcco  ship  sonars 
has  ltd  to  tho  proposed  dovclopr.icnt  of  lar£;e,  bi^h-povforcd  ^Confonoal/rianar  Ariay 
sonar  systems.  Those  systems  enhance  detection  through  increased  surveillance 
coverage,  higher  resolution  capability  and  higher  transmitted  pouer  lovely 

Much  of  the  potential  enhancement  of  detection  is  duo  to  the  in.^re.'ised 
volume  of  data  inherent  to  the  system.  Thus  parallel  to  the  physical  development 
of  an  array,  means  must  be  devised  to  process  and  extract  target  information 
from  very  large  volumes  of  data. 

Vli.'o  distinct  aspects  of  data  processing  are  evident.  First,  reduction  of 
data  volume  vilhin  a single  sonar  ping  interval  and  second,  the  intiroping  processing 
of  data  collected  over  several  ping  intervals.  The  multiping  processing  may  bo 
consklorcd  as  the  digital  automation  of  tho  detection  function  normally  carried 
out  by  a sorisr  opox'ator.  It  involves  a simple  track  and  implies  a degree  of 
classification  but  is  to  be  distinguished  from  tho  accurate  ti’ack  and  detailed 
classification  f unctions. 

The  current  ropoi't  is  tho  third  in  a series  on  multi-ping  processing.  The 
first  (1)  considered  an  H/U  processor;  one  which  required  M geometrically  correlated 
single  ping  threshold  crossings  vxithin  a span  of  N ping  intervals.  The  performance 
of  this  processor'  was  below  expectations  and  was  adjudged  inadequate  due  to  the 
poor  use  of  available  amplitude  information.  In  tho  second  report  (2),  a foi-m- 
ulation  of  an  optimum  multi-ping  detector  is  undertaken  which  close.ly  follo;:s 
veil  dovelopovd  radar  design  techniques  for  the  processing  of  pulse  trains.  The 
formulation  takes  into  account  tho  stationarity  diffore.accs  between  tho  sonar 
problem  and  tho  ideal  radar  problem. 

The  intent  of  the  current  report  v.ms  to  extend  tho  work  of  the  second  to 
include  a qiiantitativo  porforinanco  evaluation  using  current  best  estimated  of 
tho  statistical  distribution  of  sonar  signals.  Signal  models  and  their  associated 
performanco  equations  wore  developed  but  the  effort  was  tormi.natcd  prior  to 
programming  of  the  equations  on  a digital  computer.  This  repoi*t  records  tho  work 
accomplished  and  serves  as  a possible  starting  point  for  future  investigations. 


T 
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II.  TH!i  OITIHUM  I‘liOCtr^t:OR 

A dotailrd  doiivat j o;i  of  the  form  of  tho  theoretical  optiniiini  data  procosf.or 
for  a sonar  system  v.as  derived  in  (2).  Portions  of  that  work  eivc  repeated  here 
end  may  be  recognized  as  the  classical  devolopaicnt  foi'  testing  statistical 
hypotheses  given  K samples  from  a unive>rso  with  an  ass’.imed  distribution  but 
unknov.'n  paramejtcrs . 

Lot  D be  tho  total  numbe;r  of  discrete  locations  under  survoillanco.  Lot  A 
be  tho  number  of  locations  acci'ssiblo  to  the  target  botv.’oon  successive  ping 
Intervals.  Then  tho  number  of  possible  target  vmtlis  after  tv:o  intervals  is  DV. 
Following  tho  pingi  there  are  DA  possible  pat’ns  v.hich  tho  target  may  have 
traveled  through  the  ^urvoillanco  area.  Tho  data  processor  is  to  dccido  between 
tho  altcrjratc  hypotheses: 

H(l)  : target  is  present  v;ith  path  1 

H(DA^  ! target  is  present  with  path 

H(DA^’”^+l):  target  is  absent. 

The  data  from  ping  k may  bo  summarized  in  voctor  form  as 

• • • • • ^ k“  1|  2,  ....  N 

whore  is  tho  single  ping  output  corrosponding  to  tho  location  on  ping  k. 

It  is  desired  to  process  tho  data  from  successive  pdngs  in  a way  which 
maximizes  the  probability  of  detection  for  a specified  average  number  of  false 
alarms.  It  is  well  knov/n  that  tho  data  processor  which  accomplishes  this  objective 
oust  compute  the  set  of  a posteriori  probability  density  functions 

L(Z,N,j)  = p(Zj Zj^l  H(j))  l.C.^1,  ... 

and  compare  those  with  a threshold  C which  is  inversely  weighted  by  the  probability 
that  Il(j)  is  true.  That  is,  if 

L(Z.N.j)  > C_ 


(1) 
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nccopL  H(j),  olhorwino  H(j)  is  rojcoicd.  Tho  sinclo  pin^  processor  ouLpuls 
corresponding  to  different  locations  and/or  different  pin^s  arc  assumed  to  bo 
independent  allowing;  the  representation 


N D • 

L(z.N.j)  = 7r  rr  p(\i  | h(j)) 

k=l  i=l 


(2) 


The  single  ping  processor  consists  of  a bank  of  D matched  f ilter/cnvclcpc 
detector  combinations,  one  for  each  location  in  the  survoiliancc  area.  V/ith 
target  absent,  the  output  is  assumed  to  have  a Rayleigh  distribution  given  in 
tho  usual  form  as 


p(R)  = R_  exp 

; 


(3) 


whore 


R = amplitude  of  tho  voltage  onvelopo 

~ mean  square  value  of  noise  voltage-  VJith  target  present,  tho  output  is 
assumed  to  have  a Rico  distribution  (non-scintillating  target  model)  and  is  gii 
by 

P(K)  = _R_  exp  - 1 ( _RP_  ) R i 0 ( 

• I 2 Y ■ J ^ l|/„ 


where 

R = amplitude  of  tlio  sine  wave  signal 

I = modified  Bessel  function  of  tho  first  kind 

o 

The  optimum  processor  must  compute  tho  expression  given  by  equation  ? whore 
P(\i  I H(J))  is  Rice  distributed  if  the  indices  k and  i correspond  to  the  true 
target  location  under  hypotnesis  H(j),  and  is  Rayleigh  distributed  otherwise. 
Substituting  equations  3 nnd  h into  2 yields  for  the  first  DA  hypotheses 

rr  I TT  Sfi  r,  [.  fLiy.Lt.pi]  i/fyusii' 

k--ij  1=1  Yo  I 2 „i/  i J ^ ^''o  ' 

. i 5^  a(kj)  > 


v.'horci  a(k,j)  is  tko  tiuo  locaLion  of  the  targot  on  pin,:^  k ur.dcr  b;,polhcsis  j. 
Iho  tost,  function  t.iay  bo  rowi’iCton  tiioro  sir;ply  as 

= e>^i>  [-  ^7  ] f f 1 > 

L J ({/„  =>11 J ° I v/„  / 


i,/. 


Those  factors  of  L(R,N,j)  v;hich  ai-o  the  sano  for  each  hypothesis  need  not.  be 
conputcel  since  they  x;ill  haa’’o  no  effect  on  the  relative  p;a£;nitudo  of  L(R,N,j). 
Simiiarly,  tho  constant  factor  exp  KP^/k  r:ay  bo  included  in  the  cor;pari; 

circuitry.  Thus  tho  tost  reduces  to: 

Accept  H(j)  if 


li'n  / 


> 


If  tho  test  function  does  not  exceed  the  threshold  for  any  possible  path,  the.i 
tho  target  is  declared  absent. 

Kov;  that  tho  optiral  tost  function  has  been  obtained,  any  n;onotonic  function 
of  it  nay  be  used  just  as  vroll.  It  is  convenient  to  consider  the  logaiithin  of 
h(K,h,j)  vhich  yields  tho  form 


I 

Ir.  L(K,K,j)  + In  > In  C 


If  we  assume  all  paths  equally  likely,  then  tho  path  probability  contributes 
equally  to  all  tests  and  nay  bo  ignored.  Thct  final  decision  criteria  becomes : 
Accept  H(j)  if 


Vo  I 


> In  C 


Thus  tho  optimum  processor  takes  the  sum  of  R variates  where  each  variate  has 
boon  processed  by  the  function  Indicatt-xl  in  equation  9-  Tor  siiall  signal  to 
noise  ratio,  the  following  expansions  may  bo  applied 

I (X)  = 1 4 i 
o 


In  (1  + i X^)  i X^ 


_ -a 
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Insortiti"  thor.c  oxjiansiojis  in  tlic  tc-rjt  yinl/lc  t’no  losult 

\ J.2  > constant  (ll) 

16l|//  L "•«(«) 


that  is,  tho  optinal  pi-ocossor  cr.n  bo  roplacc-d  by  a,  square  law  detector  with  a 
linear  integrator  for  small  signals. 


_6.. 


m i'w(rOi{i-:\;;cE  mi.UATic:;  - stevuy  sig.\’als 

Tho  oj'tirnun  procosror  lia.s  been  chcv.’n  to  consist  of  tho  linear  integration 
of  sono  fonclion  of  tho  onvolopo  of  tho  single-  ping  pi'ocos.sor  output.  It  is 
conv'c-niont  to  take  tho  output  as  that  of  a square  lav:  detector  and  normalize  with 
respect  to  tho  moan  square  noise  voltage  tov.'ard  this  end,  let 


y = 

(1?) 

X ='  = S 

^Vo  ® 


The  quantity  x may  be  identified  i.'ith  tho  signal-to-noise  pov.'cr  ratio  commonly 
used  in  detection  analyses.  In  terms  of  tho  now  variables,  the  single  ping 
output  distribution  of  intcrost  are: 

For  noise  alone 

=-■  exp  yj  y > 0 (I3) 


For  signal  plus  noiso  — steady  signal 

fjCy)  = exp  [-y-xj  Ijj  (2'V^  ) y>0 


(1^) 


lor  a single  ping,  tho  probability  of  false  alarm  and  detection  are  equal  to  the 
probability  of  exceeding  a preset  threshold  with  the  target  absent  and  present 
respectively.  In  equation  form,  tho  probabilities  are; 

o/> 

ti  - j Oy 

c 


’’d  - J hW  S 

e 


(15) 


For  the  sum  of  N variables,  lot 

N 

L y* 

1=1  ^ 


y = 


(16) 


The  distribuLion  functions  of  Iviualion  15  roplaced  by  the  dir, lributi.cn  function 
for  Y.  Thir  is  nost  easily  obtnintd  through  tho  uso  of  characteristic  function.'.: 
vhorein  tho  characteristic  function  of  tho  su;i  of  li  variates  is  equal  to  the 
product  of  tl)oir  individual  cliaractoi'istic  functioiis. 

For  signal  plus  noise,  tho  c>!aractoristic  function  of  Dquation  ih  is  given  by 


CjCp)  j expj-x-yj  ) oxp  Jpyj 


T)io  integral  my  bo  obtained  fron:  pair  of  Canpbcll  and  Foster  [ 3] . The 

integrals  of  this  leforcnco  use  oxp  -pyjfor  the  first  integration  and  exp 
foi-  tho  inverse  transforn.  This  Siiriio  practioo  is  used  boro  to  al.l.ov,’  direct  use 
of  the  tables  of  reference  3*  Inserting  the  -p  in  iquation  I7  yields  the  function 

Cj(p)  = exp  l^-xj  exp  |^x/(p^l)j  (18) 

Tho  characteristic  function  for  the  sura  of  M variates  is  then 

C:j^(p)  = [ (19) 

“ 1 oxp  ®^P  ^Nx/(p'il)j 

(p+1)^ 

From  pair  65O.O  ^3j  # the  donsity  distribution  is  seen  to  be 

Nd. 

f^(T)  -■=  ^ oxp[-y-Nx]  I^,_j  (2-yiTxY)  (20) 

Tho  probability  of  detection  is  given  by 

VN)  = J fj^  (y)  dY  (21) 

c 

This  intogiTtl  is  not  soluablo  in  terms  of  well  known  functions.  It  is  a siecial 
case  of  the  inconiploto  Toronto  function  which  is  dorined  as 


Ty(m,n,r)  = 2r 


n-m+1 


I (2rt)  dt 

n 
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In  trims  ( 

af  tliis  fund 

:.ion,  tlie 

dr  Led  ion  pr 

ob'ibility  is 

yn)  - 

1 - T 
Vc 

(2:J-1,  n-1. 

V^Nx  ) 

(23) 

Curves  of 

tho  incomplete  Toi'on 

to  function  a 

ro  available 

in  refei'oncc 

and 

numerical 

means  of  evaluating 

the  function 

are  available 

in  the  HMi® 

piogramming 

library. 

For  noise  al 

Lone,  an 

identical  procedure  may  be 

carried  out 

vrith  tho 

re.-ulUant 

equations : 

Roj)cating 

Equation  13 

fjy)  ^ 

^ £->:p 

li'om  Kefe: 

I'cncQ  3 pair 

Ki'i  or  d 

irect  iniegra 

tic>n 

Ci(p)  ^ 

1 

(24) 

1-tp 

Then 

Cj.(p)  = 

(l-'p)*'^ 

(25) 

From  Reference  3 pcir 

^!31 

r"-»  ...  ! J 

f (Y)  = 

(N- 

1)1 

(2.6) 

» 


♦ 


Tho  probr.bility  of  false  alarm 


p^OO  = 


f f (Y)  dY  = 1 - I 

f C , N-1 

C 

Ivn 

(2?) 


vhorc  I (. ) is  Pearsons 
by  tho  expression 

I(u,p)  = 


form  of  tho  incomplete  Gamma  Function  designed 


o 


(23) 


Tho  pi-obability  of  false  alarm  as  given  by  F/]uation  27  pertains  to  a particular 

N-1 

possible  target  path.  Tho  number  of  such  pc»ths  was  defin'd  in  Section  II  as  DA  “ . 
For  largo  A,  the  processor  test  functions,  L(R,N,j),  are  approxii..at  ely  independent 
f~  of  cacli  other  and  tho  probability  that  n of  those  result  in  a false  alarm  is 


t 
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{”rA(!;)  - = 


r 1*" 

r 1 

1 1 - rpii): 

(7-9) 

\ n / I 

^ j 

L J 

Thf'  .ivci'ar^o  nuiibc  r of  fal:;c  alarms  is 


* r -\ 

n \\  <;FA(:J)  = nj> 


n=l 


(30) 


The  equation  for  lA(r;J  is  roen  to  bo  equivalent  to  the  oqirossion  for  th; 
iscan  of  a binomial  d istributioii  v;hich  has  the  closed  foim  lesult 


FA(N 


DA^*"'  P,(ti) 
A 


(31) 


V.'ith  this  definition  of  average  mimbor  of  false  alarr:s,  ’.:o  note  that 

.N-1 

PA  appears  as  a simple  multiplier.  If  v.’o  construct  a plot  of  — or 
alternatively  signal-to-noisc  ratio  requirc-d  for  some  fixe^d  value  of  P„  -- 
versus  PA(iry,  vo  need  only  plot  against  P vrith  DA  ” as  a scale  factor. 

A 

liquations  23  and  7:j  were  evaluatrd  by  oxcercising  existing  in-house  sui)- 
routinos  for  tho  indicated  functions  for  soloctc-d  throshold  values  C end  S/ii 
ratio  X.  Curves  of  Pj-j  voi’sus  P^  for  fixed  values  of  I.’  and  x,  plotted  for 
common  values  of  C,  are  given  in  Figures  1 thru  6.  A cross  plot  of  this  data 
is  given  in  Figure  ? wliich  sho'xs  the  S/N  ratio  in  db  required  to  achieve  a 
detection  probability  of  0.9  as  a function  of  number  of  pings  processed  and 
probability  o^  false  alarm. 
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]V  COMPAIiKO:!  Or  OiTi:-DI-!  AI.'D  M/U  rKCCfi;.SOK--nO:i-SCl;:TIMATL’G  TyUtGLTf) 

A dcilnilcJ  troatniont  of  an  M/l.'  proconror  is  given  in  Koferenco  1.  This 
processor  operates  under  a target  dotoction  criteria  in  wliich  M goor.otrically 
correlated  single  ping  threshold  crossings  are  required  vithin  a sjvin  of  N 
ping  i ntc;rvals.  A simpler  treatment  of  that  problcin  is  given  hero  to  allov; 
comparison  of  the  K/K  processor  with  tiio  optimim  pj’occssor. 

Tho  cumulativo  multi-ping  probability  of  detection  and  false  alarm  are 
related  to  the  single  ping  values  by  tlic  oqaations  : 


Pu(N)  ••= 


£ I : I I •-  '“'”1' 


/N\  1,  , 

L k 


As  in  the  optimum  processor,  the  average  number  of  false  alarms  is  given 


FACO  = P^(N) 


In  tho  comparison  -which  follows,  N is  fixed  at  five  ping  intervals  and 
the  value  of  D and_A  apply  equally  well  to  both  optimum  and  K/K  processors. 

Thus,  a full  comparison  of  tho  tx'.-o  types  of  processors  may  be  accomplished 
without  considering  the  ©xplicite  values  of  D and  A. 

The  single  ping  probability  of  dotoction  and  false  alarm  are  related  to  each 
other  by  their  common  threshold.  From  Equation  l^A,  v;o  write 


ppi) 


■ r 


thus  the  thi-cshold  may  be  written  directly  in  terms  of  false  probability  a: 


c = In  II 


Applying  this  thioshold  to  Equation  yields  for  probability  of  detection; 


p„(o  = 


(P^GJ)  dy 


'(P.n^) 


o- 


This  intf-^ral  has  no  ’losod  form  arid  is  pcnr: rally  cy.prtir.vi'.'i  as  Karciins  Q 
function  { \ 


i/r/i) 


The  computational  proccduro  for  tho  coi'-piriscn  included  hero  is  as  follovrs: 
1)  Set  N equal  to  5 find  Pjj(5)  equal  to  0.9. 

?)  Usin<;  K-qintion  J7.  find  equal  to  1 thi-ou^h  5« 

3)  for  each  M,  use  Equation  33  to  find  the  S/li  i-atio  required  for  selected 
values  of  T (1  ) 

/V 

4)  Using  Equation  33  find  tho  P^(5)’s  v.’hich  correspond  to  the  of 

step  h. 

5)  Wot  the  S/li  required  versus  P.  (5)  for  each  M 

’ A 

The  results  of  these  calculations  together  with  tho  optimum  processor  for 
N equal  to  five  arc  given  in  Piguro  8.  In  general,  the  optimum  shov/s  about  a 
two  db  advanUigo  over  tho  M/K  processor. 


* 
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SCKiTfLUA.TJIlG  TA!;GLT 


Tho  nnaly.'jis  to  thii;  point  utilizer,  n fixfj  S/N  ratio  at  tiio  input  to  the 
rocoiver.  That  thi.a  ir  not  t’no  ease  is  woll  knov.-n  and  suitable  no:Iif ioations 
must  be  included  in  the  analysis  before  it  can  bo  con';ido''cd  realistic. 

t 

Tho  most  accepteJ  moriol  roprosontative  of  variation  in  si(;nal-to-noiso 
Is  tho  loc-nornal  di.stribution.  This  model  has  been  oxtonsively  used  in 
si£;nal  pin^  analyses  and  does  have  justification  from  tot'n  theoretical  and 
expcrimnntal  standpoints. 

Accepting  tho  lop,-no?'mal  distribution  as  tho  best  tarj^ot  scintillation 
model  .allov/s  v-’niVtin^  the  distribution  in  tho  form: 


f(Z)  “ ^ exp 

vhoro  Z = S/i'i'  povrer  ratio  in  db 


2^Z 


(39) 


Z = moan  of  Z 

Z - stand.ard  deviation  of  Z 

The  general  formulation  of  porformanco  analysis  requires  expr-esfing  the 
distribution  in  terns  of  .S/N  povrer  ratio  in  volts.  The  change  of  variable  required 
is : 


Z = 10  log  X - 10  log  e In  X 


(^0) 


Completing  tho  formal  change  of  variable  yields 


f(x)  = 


V27r<^': 


X X 


whore 


X = 


10  log  o 


(^1) 


Attempts  at  forcing  tho  distribution  of  E:iuation  41  into  the  analysis  of 

Section  III  proved  hopeless.  Tho  principle  difficulty  v:as  tho  appearance  of 
2 

(inx)  terms  wliich  cannot  bo  conveniently  represented  in  series  form. 

To  circumvent  tho  difficulty,  a distribution  v;as  assumed  vdrich  had  character- 
istics similar  to  tho  log  normal  and  the  analysis  completed  using  the  assumed 
distribution. 


r 
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A fnnily  of  ccititill atin^  l^r^ot  nc.-.lcl;; 
vith  the  {jotioivil  equation. 


pror.cntrd  by  Sv.’crlin^ 


w (x  , X ) 


1 K 

Tk-07  = 


exp 

r-  Kx] 

1 J 1 

1 

1 _ 

!.  X . 

('12) 


Tho  equation  is  a Chi  square  (iistribulion  vith  ?K  do2roos  of  froc.doio  vrhcrc  x 
and  X have  tho  samo  intorprotation  as  in  liquation  41. 

Tho  best  value  of  K proved  to  be  3/2  — a Chi  sqiiaic  distribution  with 
3 degrees  of  freedom.  For  thbvaluo  of  K,  tho  equation  for  S/'.l  distribution 
bocoiiios  ' 


v(x,x)  “ , ? 

V/r 


1 

X ^ exp 


I 

2x  J 


(43) 


A comp.aritivo  plot  of  Equations  4l  and  43  aro  shown  in  Ticui’c  o.  Ibc  uso 
of  Equation  43  as  tho  S/iJ  distribution  appeal,  rocssonable  particulai’ly  in  view 
of  tho  uncertainties  inherent  to  Equation  41. 
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Tlio  f. plurj  noi.'ir  cl i retribution  for  fixed  t, ijjiVii -to-noir;c  ratio  wan 
given  by  Ejuation  an 

f(y  1 x)  exp  ^-y-xj  (?  ^ y.y  ) 

Tho  a.nnumrd  cl intribution  for  signal-to-noinc  variability  vran  given  by  Equation 


43  an 


v(x)  - J> 

iir 


a_ 

? X 


3/2 


exp 


- 

2x 


Given  thene  clinti  ihutionn , the  unconditi  onal  probabii.ity  of  y in  general 
form  in  given  by' 

f(y)  ~ J f(y|  jc)  (4-0 


For  our  problem,  tho  rcnultant  equation  is 

f(y)  =-  b cxp(«y)  J x^  exp  (-ax)  I (2  "^/xy  ) dy 


(43) 


whore 


= _3_ 

X 


7k)  ‘ 


a = 1 + 

2x 

The  Bosr>ol  function  may  bo  written  in  serien  form  as 

ei‘> 

(2  Ay'  ) = L 

k"0 

which  allows  ro-writtin-^  tho  integral  of  K'liMtion  as 

r 

JM.'!  / X ''■t 

^2 


(klf 


Y. 

k==o 


(kO' 


exp 


i-I 


dx 


(46) 


(4?) 


To  evaluate  tiio  integral  of  fiquaiion  47i  we  note  tho  Gamma  Function  in  defined 

t>y 

d t = f'(Z)  (4S) 


,Z-1  -t 


<1 


~?>i- 


K-ikin"  fiti  appj'onri.-'.tc  cli  of  variable-,  the  integral  of  B-^uatlori  hj  n.ay  be 

written  in  the  clo55cd  form 

k»'i/? 

[--—'j  ' f’(k  ^ 3/??) 

The  i-ncotvUtional  dir.tribution  of  y thus  hoco!.!05; 


f(y)  =--  b exp  (-y) 


Fejualio.i  U')  pc-rtainr  to  the  distribution  of  a single  pulse.  The  distribution 
for  the  sen  of  b pulses  is  obtained  in  a rxinner  analogous  to  that  used  in 
Section  in. 

To  siriplify  notation  let 


k kl 


k ^ 3/2 


Equatio.n  r.ay  non  be  v.'i  ittc.n  as 

f(y)  = exp  f-y]  (51) 

k--o 

I'rotn  Campbell  and  Forster,  pair  h31,  the  characteristic  function  for  one  pulse 
is 


c,(p)  - 

' to  '>'*»> 

For  the  sum  of  I>  pulses,  we  write 


c.,(p)  = ) ’ • Y ^ 

r*  1^1  4 u 


kj-o  k^-o 


J--J  . .V  kl  + k._,  ...  k.  + 

k,.-o  (p*i;  ^ N 


Urin^^  p:iir  ^^3^  inverno  trancrorin  yiolds  tho  dir^t-ribution  of  the  suiu  as 


f(Y)  = 


\ ' . . . . y ^kl  •••  •^kh exp  [-Tf]  "y  ^ 

/ ->  I ->  41,  U 4 K _ -1  1 I 


ki-'k2...k.^-iM-l 


I V ^ (k.^k.,  ...  k^  N - 1)1 


Tho  probability  of  dotcctLor.  ir  obtn in- 1 by  integrating  Piqnation  5^  froi.i 
Uk!  iliroijhold  C to  infinity.  Taking  Uif  Y torn3  only  wo  havo 


C«’ 

J oxp  [ -y| 

rl' 


..vL  • dy 


(:i5) 


vhero 

r ^ kj  + ^ -•  N - 1 

This  intf-gral  has  the  s.anio  forfl  as  Poai'son  Inco;nploto  Gapirn  I’onotion  •.rhich 
is  defined  by  liquation  ?3.  Symbolically  tbo  integral  is  rc;p’’csentf  1 by 


1 - I 

IVi-n 

Tho  final  losult  is 


. r 


pj^oo  - y y y ( •••  \ij  ) 

kj=o  k2~o  kjj“o 

vhero 

r = k,  •*  k,j  . . . k.  ^ K - 1 
12  N 

T(  • ) ~ Pearson  Incomplete  Gamma  Function 

A.  = b r^^3/p.) 


ki 


b = ^ 

X 

1 + ^ 

2x 


k.  + 3/2 

a 1 
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zU 


1 - I 


■1 


(56) 


Tho  false  alarm  conputation  is  identical  to  that  given  in  Section  III. 

A comparison  of  optimum  and  M/K  processors  may  bo  made  in  a manner  analogous 
to  that  of  Section  IV.  TTie  only  change  required  is  a change  in  the  signal  to 
noise  distribution  function  given  by  Equation  51* 

Piaking  this  change,  Equation  37  Section  IV  bocemes 


O^' 


Ih-  iiitfjril  is  ar;ain  r«  co-nir-ibln  ns  Pc.u.'.ons  IiiCo;r.j,lrfio  Ga'-irsi  I'ur.cLion 
vith  Uiu  fiti-i]  result. 
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Uso  of  on  cptinun  proccsror  ir.:plic";  tho  prococrin^^  of  cnormour.  amounts  of 
dain.  It  is  ai)parcnt  that  a practita]  processor  nust  contain  data  roJucing 
tcchnlqurs  to  allovj  coraputational*  feasibility.  Ono  such  technique  is  t'nrosholdin^; 
the  single  pin?;  output  prior  to  multi-pin^  processing.  Tho  njulti-ping  pjocessor 
then  concerns  itself  only  with  outputs  v/hich  ozcotd  a pi-osct  volta^io  level. 

Tho  overall  process  bocomor.  a double  threshold  device  and  perfor.'^nncc  is  a 
function  of  tho  setting  of  the  tv;o  thresholds. 

The  effect  of  single  ping  thresholding  on  tho  detector  output  is  piciorially 
shown  bol  ov; 


Ibat  is,  throsholdini;  places  an  impulse  at  y equal  to  zoro^cqual  to  the 
probability  that  y is  less  than  C^. 

Tho  distribution  function  may  be  written  as 

c(y)  = ^ “ (y-Cj) 


(57) 


vhcro 


£ (y)  ~ Dirac  Delta  Function 


u(y-Cj)=  I 


0 

1 


y^  Cj 


y ^ c, 


For  noise  alone,  tho  distribution  was  p"cvlously  shown  to  be  given  by 
f(y)  - exp  ^-yj 
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lov  a fixed  valuo  of  C^.  ^ho  va}i;o  of  ^ bocoir.er  a constant 


of  the  prolilon  and  i.s  equal  to 

Cl 


I 


exp  [^-y]  dy 


(y3) 


= 1 - oxp  J-Cjj 

The  threshold  distribution  function  nov;  bocories 
e(y)  = Aj  C(y)  + exp  ^-yj  u (y-C^ ) 


(5?) 


The  distribution  for  the  sum  of  N pulses  is  a^^in  derived  usin^  characteristic 
functions.  The  Fourier  transform  of  tho  delta  function  is  unity.  Tho  transforra 
of  tho  remaining  term  is  obtained  by  direct  integration  as  follovrs 


J exp  J-y]  u(y-C^)  oxp  ^-pyj  dy 
o 

= J exp  ^-(1  + p)  yj  dy 
c 

\ 

®^p  [-  (i+p)  cj 

Tho  characteristic  function  for  a singic  pulso  may  be  written  as 
C.  (p)  = A.  1 exp  f-  (l-*p)  cj 


(60) 


(61) 


Tho  charactoristic  function  for  tho  sum  of  F pulses  is  tho  N-fold  product  of 
Equation  6l.  Using  tho  binomial  thoorom,  tho  resultant  nay  bo  written  in  finite 
sories  form  as 

F 


**  nr 

L 0 


N-k 


‘1 


i exp  [ -(l^p)  k C 1 iC?.) 

^1^7)^  L IJ 
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Tho  cxjionciiti-il.  in  p is  simply  a cl  1 spl •loomonl-  — i.o.  from  Campboll  and 
Forstor,  Pair  207 

oxp  f-  pcoj  F ->  G (g  - Cjj) 

Tlio  invci'so  transfoim  is  a£^;ain  obtained  from  Cinpboll  and  iorstc-r,  pair 
v;ith  thn  rosult. 


k--o 


(y-kCj 

“TCiT” 


(63) 


Tho  probability  of  falso  alarm  is  obtained  from  Equation  63  by  intfcratine; 
from  tho  multi-pina  threshold  to  infinity.  For  tho  terms  containin3  Y,  tho 
integral  is 


(Y  - kCj oxp  [ -Y]  dY 

"Tk“)r 


(6^0 


To  perform  tho  integration,  first  make  tho  change  of  variable 
, Z = y - kc 

vhich  results  in  tho  expression 


oxp 


C„-kCi 


yk-1 

~(k-T)r 


oxp 


f-Zl  d! 


(65) 


Tho  latter  integral  has  the  form  of  Pearsons  Incomplete  Gamma  Function  defined 
by  liquation  28  yielding  for  the  integration 


oxp 


, k 


Vit 


-11 


Tho  final  expression  for  probability  of  falso  alarm  is 
I N I N_v 

[-xh] 


V t N \ N-k 

IJ  *. 


oxp 


1 - I , k-1 


/k 


vdicro 


C,  - kC, 


Aj  - 1 - exp  [-C;^ 


(66) 
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CoPiput.,nli on  of  the  averaco  minbcr  of  falr-o  alainu  ir.  identical  to  tho 
earlier  dovolop?;ont  in  Section  III  and  i':  given  by  Equation  31* 

When  ta  steady  target  is  present,  the  distriiiution  of  y is  given  by 

f(y)  ~ exp  1^  (?-  ) 

The  probability  that  y is  less  than  tho  ninglo  ping  threshold  value  is  easily 
shown  to  be  exprossiblo  in  terms  of  Karcuns  Q function 

?r  <{f(y)  = ^2  = 1 - Q(b,^.  ‘/iCj  ) 

Tho  tiiresholded  distriVmtion  then  has  tho  form 

c(y)  " ^2  ^ uCy-'^j)  (6?) 


bsing  the  series  expansion  for  I given  by  E-quation  yields 


e(y)  ^2  aCv)  ^ ^ 


hx)- 


exp 


k=o 


(kr)^ 


(«i) 


A general  dovolopncnt  for  computing  tho  probability  of  detection  for 
the  sura  of  " pulses  given  a single  ping  distribution  in  tho  form  of  Equation  C8 
is  given  in  Appendix  A.  Direct  appliration  of  that  work  yields  tho  result. 
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Foi-  .1  rciriti]  tnrr;ot.,  Uio  distribution  as  nivtii  by  filiation  51  is 


f(y)  “ ) „k  y*'  oxp  [ ~y] 

U >^1  ^ J 


Tho  probability  that  f(y)  is  loss  tlian  tbo  cinglo  ping  Uirtshold  is 

c,  k r 1 

V t-yn  l_v' 


- A, 


^ C k 

X \ r - »>• 

k=o 


\ ' I -;b . k 

Vk^l 


(70) 


vhoro  !(*)  is  the  Inconplote  Gamma  Function  of  Equation  28.  Tho  thresholdcd 
function  thus  becomes 


g(y)  = A^  C(y)  )2  ^ y*'  exp  [-y]  u(y-C^) 


k=o 


(71) 


Applying  tho  general  development  of  Appendix  / tc  Equation  71  yields 
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For  ^icvcri'.l  assumed  i-iuinl-plus-noiac'  distribuiicnr; , the  cx^nridcd  foriri 
of  tho  distribution  has  tbo  ^^cnoral  equation 


e(y) 


t.'  • 

Aj  s(y')  ^ t*"^] 

k-6 


(A-l) 


vhoro  y = normalized  output  of  a square  law  dotootor 
(y)~  unit  impulse  at  y ~ 0 
u(y-t^)”  unit  step  at  y - c 

The  probability  of  detection  for  the  sun  of  N pulses  from  the  above  distribution 
function  is  desired.  This  is  obtained  by  the  usual  application  of  characteristic 
functions.  Ln  v:hat  follows,  the  characteristic  function  is  taken  as  the  Fourier 
transform  allow'ing  use  of  Campbell  and  Forster  Fourier  Transforoi  Tables. 

The  Fourifu'  Transform  of  the  impulse  function  is  unity.  The  transform  of 
the  second  torn  of  Fquation  A-l  procedes  as  follov.'s  : 

Considtirinc  the  y terms  v;ithin  the  summation,  the  transform  in  integral  form 


is 


■ jy^  exp^-yj  u (y-c)  exp  [^-pyj  dy 


(A-2) 


Tho  unit  step  at  y = c may  be  eliminated  by  sotting  the  lower  limit  at  C 


I- 


(A-3) 


y''  exp  J-yJ  exp  ^-py]  dy 

Since  Fouiier  transforms  are  available  with  lower  limit  equal  to  zero,  wo 
now  make  tho  change  of  variable 
Z = y-c 
dZ  *=  dy 

which  yields  tho  equivalent  expression 

^k 


I 


(Z  ^C^  exp  [-(Z+C)]  exp  [-p(Z^C)]  dZ  (A-4) 

V 

Tho  binomial  torm  my  now  bo  written  as  a finite  series  with  tho  result 
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(J  c z^ 

exp  ^-Zj  exp  ^-(I'p)  Cj 

o.p  [ -pZ^  dZ 
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(A-5) 
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Cami'l'‘''T]  and  dorstor,  {lair  hjl,  proaonl  tlio  transfon 
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E > 0 


(p+  y 


Taking',  the  inlc£;ral  inside  the  sif.nr.:'iti on  and  ajolyintj  pair  ^01  yields 


exp  [-Z]  ^ 


The  P'ourlcr  trjnsforin  of  Equation  A-1  ray  novr  be  v;ritten  a' 


Cj(p)  = Aj  + y 
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Tho  Fourier  transform  for  the  sum  of  K pulses  is  the  K-fold  multiplication 
of  Equation  A-?.  Sjanbolically  treating  the  right  hand  sido  of  Equation  A-?  as 
tho  SUM  of  tv;o  terms,  wo  write 


I PI 


N-tn  , m 
a b 


Inserting  tho  appropriate  expressions  for  a and  b yields 
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(A-8) 


Ibo  next  stop  is  to  take  tho  inverse  transform.  The  only  terms  wo  need 
consider  aro  those  containing  p.  Theso  aro: 
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exD -rmc ' 
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Tlio  exponential  term  ropresoats  a dir.placcmont,  i.o. 

-pYo 

c(p)  c £(y-Yo) 


The  deno;ninator  of  Fquation  A-9  nas  tho  sanio  forr:  an  tho  trnnnfor.a  p;.ir  given 
in  E>^uatlon  A-6.  The  resultant  transfonn  of  Equation  A-9  is 


exp  [-pmc]  j (Y-r.c)^-'  exp  ^ (-Y-tnc ) 1 

— Tr-TJl  ^ 


(p+ir 

The  integral  of  g(Y)  from  a tnulti-ping  threshold,  to  infinity  yields 
the  probability  of  detection.  To  carry  out  the  integration  v;o  may  restrict 
ouj'solves  to  those  terms  containing  Y as  given  by  Equation  A-10.  The  desired 
integral  is 

I (Y-mc)'"^  exp  l-Y^  dY 


^rej  (r-l)l 


To  carry  out  tho  integration,  first  make  tho  change  of  variable 
Z = y - me 
dZ  = dy 

which  yields  tho  equivalent  intcgi’al 


f 
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Cjj-  me 


The  latter  integral  may  bo  identified  with  Pearson’s  Incomplete  Gamma  Euriction 

which  is  defined  as  

u 

I (u.p)-  f 


In  tonus  of  this  function,  Equation  A-12  becomes 


C„  - me 
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Tho  final  oxproi;£,irjn  for  probability  of  dotnction  is 
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